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^ ■ Abstract 

^n ; 

t"^ , We study the effect of squark generation mixing on gluino production and decays 

C^ I at LHC in the Minimal Supersymmetric Standard Model (MSSM) for the case that 

• ■ the gluino is lighter than all squarks and dominantly decays into three particles, 

L». , g — ;• qqXk^ QQ'Xi ■ We assume mixing between the second and the third squark 

T-H I generations in the up-type and down-type squark sectors. We show that this mixing 

^~~^ • can lead to very large branching ratios of the quark-flavour violating gluino three- 

^ . body decays despite the strong constraints on quark-flavour violation (QFV) from 

>v>( I the experimental data on B mesons. We also show that the QFV gluino decay 

;h ' branching ratios are very sensitive not only to the generation mixing in the squark 

. - - I sector, but also to the parameters of the neutralino and chargino sectors. We show 

that the branching ratio of the QFV gluino decay g — )■ ct{ct)xi can go up to ~ 40%. 

Analogously, that of the QFV decay g — )■ sb{sb)xi can reach ~ 35%. We find that 

the rates of the resulting QFV signatures, such as pp — t- ttccEJp^^, can be significant 

at LHC. This could have an important influence on the gluino searches at LHC. 



1 Introduction 

The flavour structure of tlie quark sector is very well described by the Cabibbo-Kobayashi- 
Maskawa (CKM) mixing matrix, which is the only source of quark flavour violation (QFV) 
in the Standard Model (SM). In particular, flavour changing neutral current (FCNC) 
processes, such as K'^ -^ jJ^'^ l^' , B^ — ?> fi'^fi^, B — )• Xs7, B — > Xgl'^l" etc., are strongly 
suppressed [1]. They impose strong constraints on the quark generation mixing. Any 
extension of the SM must therefore respect these constraints. 

In supersymmetric (SUSY) extensions of the SM, mixing between different quark 
flavours in the squark sector, that is not related to the CKM-matrix is also possible. 
Although the mixing between the first and second generation squarks is strongly con- 
strained, there is room for appreciable mixing between the second and third generation 
of squarks, still obeying the constraints from B meson data^. This is beyond the minimal 
flavour violation (MFV), where the only source of QFV is the mixing due to the CKM 
matrix [3, 4, 5]. 

The effects of mixing between the second and third squark generations, especially the 
mixing between top and charm squarks, have been studied in the Minimal Supersymmetric 
Standard Model (MSSM) for squark production and decays at LHC [6, 7, 8, 9, 10]. Squark 
generation mixing has been investigated in detail in the QFV decays of gluinos, g -^ 
ui,2 c ('Ui,2 c) — )■ cixi (ct Xi), where ui^2 are the lightest squark states and are mixtures of 
charm and top squarks [11]. There it is assumed that at least one squark is lighter than 
the gluino, so that the gluino decays first into a real squark (antisquark) and an antiquark 
(quark). In [11] it is shown that this leads to pronounced edge structures in the charm 
top-quark invariant mass distribution. 

In the present paper we study the QFV gluino decays in the general MSSM assuming 
that all squarks are heavier than the gluino, so that the gluino dominantly decays into 
three particles. This will give rise to a very different pattern of QFV gluino decays 
as compared to the studies in Ref. [11], due to interference effects between the various 
virtual squark exchange contributions. Moreover, the invariant mass distributions will 
have different shapes without any edge structure in contrast to the case studied in [11]. 

We study the mixing between the second and the third generations not only in the 
up-type squark sector as in Ref. [11], but also in the down- type sector. We investigate 
QFV gluino decays including those into down-type quark pair plus neutralino, such as 
g — )■ sbxi- Furthermore, we also investigate in detail the dependence of the QFV gluino 
decay branching ratios on the neutralino/chargino parameters. We take into account all 
relevant experimental constraints on the MSSM parameters from B physics and searches 
for Higgs bosons and SUSY particles, and the theoretical constraints on the trilinear 
couplings from the vacuum stability conditions. 

Recently ATLAS and CMS performed searches for SUSY at LHC with ^/s = 7 TeV 
on the basis of total integrated luminosities of 35 pb~^ [12, 13] and ~ 1 fb~^ [14, 15]. 
They found no excess of events over the SM expectations and set limits on the squark and 



^ There could also be mixing between the right up-squark and the left top-squark which is hardly 
constrained [2]. However, we do not consider this mixing here. 



gluino masses at 95% confidence level (CL). In the simplified SUSY model with m^o ^ 
(where all SUSY particles other than the gluino and squarks of the first two generations, 
and Xi s-re decoupled by being given very heavy masses) the limit on the gluino mass 
is rrig > 800 GeV for large rrig and that on the squark mass is rUq > 850 GeV for large 
rrig [14] (see [15] also). Here, rrig is the degenerate mass of the squarks of the first two 
generations. In the context of the constrained MSSM (CMSSM) (or mSUGRA) the lower 
limit on the gluino mass is smaller than ~650 GeV for rriq > 1.5 TeV and the limit on the 
squark mass is rriq > 1.1 TeV for any rrig [15, 16] (see [14] also). Here again rriq is defined 
to be the degenerate mass of the squarks of the first two generations and the masses of 
the third generation squarks are (significantly) smaller than this rrig in this framework of 
the CMSSM (mSUGRA). Therefore we will assume a gluino mass of about 1 TeV in our 
analysis respecting these gluino mass limits. 

2 Squark mixing with flavour violation 

In the MSSM the most general form of the squark mass matrices in the super-CKM basis 
of «7o7 = (^iL, g2L, qsL, Qm, g2R, ^sr), 7 = 1, •••6, where (gi, ga, Qs) = («, c, t), (d, s, b) is [17] 

for q = u,d, where the 3x3 matrices read 

^III = M'q + Di,LL^ + ml Ml^L = VckmM^qV^^^ + D^^llI + ml 

^Irr = Ml + ^J,i?i?l + ^l Ml^^ = MIj + D^^nnl + ml (2) 

Here Mq ^^ are the hermitian soft-SUSY-breaking mass matrices of the squarks and rfiu^ 
are the diagonal mass matrices of up- and down-type quarks. Dg^^ = cos2/3m|(T3 — 
Cg siv? 6w) and Dg^nR = ^q sin^ 9w cos 2/3m|, where T§ and e^ are the isospin and electric 
charge of the quarks (squarks), respectively, and 9w is the weak mixing angle. The left- 
left blocks of the up-type and down-type squark mass matrices are related by the CKM 
matrix Vckm due to the SU(2)l symmetry. Note that VqkmMqVq^^^ ~ Mq as Vckm — 1- 
The off-diagonal blocks of eq. (1) read 



f2 _ A^2t _ ^mT 

V2 



^Irl = ^Ilr = ^T^-/i*rn.tan/3, 



MlnL = Ml^LR = ^T^-/i*m.cot/3. (3) 

Here Tjj j^ are transposes of Tj/d which are the soft-SUSY-breaking trilinear coupling 
matrices of the up-, down-type squarks: Cint 3 —{TuafsURi^ULaH^ +Tj:)ai3dRpdLaH'l), /i is 
the higgsino mass parameter, and tan/3 = v^jvx^ where v\^2 = \/2(-f^i2) are the vacuum 



Table 1: Constraints on the MSSM parameters from the B-physics experiments relevant 
mainly for the mixing between the second and the third generations of squarks, and 
from the Higgs sector. The fourth column shows constraints at 95% CL obtained by 
combining the experimental error quadratically with the theoretical uncertainty, except 
for B{Bs — > fi^fi~) and niho, which is the mass of the lighter CP-even neutral Higgs 
boson. -R|Yf7t,= BtB+^Zr+Jy^^ ~ [1 — ((mB+ tan/3)/m/^+)^]^, where mH+ is the charged 
Higgs boson mass [19]. 



Observable 


Exp. data 


Theor. uncertainty 


Constr. (95%CL) 


AMb, [ps-^] 


17.77 ±0.12 (68% CL) [1] 


±3.3 (95% CL) [20] 


17.77 ±3.31 


10^xB(&^s7) 


3.55 ± 0.26 (68% CL) [21] 


±0.23 (68% CL) [22] 


3.55 ±0.68 


10^xB{b^sl+l^) 








(/ = e or fi) 


1.60 ±0.50 (68% CL) [23] 


±0.11 (68% CL) [24] 


1.60 ±1.00 


W^xB{Bs^ fi-^fj.-) 


< 4.3 (95% CL) [25] 




<4.3 


10^xB(S+^T+z>) 


1.68 ± 0.31 (68% CL) [21] 


±0.25 (68% CL) [21] 


1.40 ±0.76 


mf,Q [GeV] 


> 114.4 (95% CL) [1, 26] 


±4.0 [27] 


> 110.4 



expectation values of the neutral Higgs fields. The squark mass matrices are diagonalized 
by the 6x6 unitary matrices R^, q = u,d, such that 



diag(m|,...,m?J 



with 



nig^ < ... < m^g 



(4) 



The physical mass eigenstates qi,i = 1, ..., 6 are given by (ji = Rf^qoa. 

In accordance with [18] we define the QFV parameters 6^^, S^^^ and 6^^^ {a ^ (3) as 
follows: 



XLL 



Ml^p/^fM^ZMl^p , 



(5) 



XuRR 



MLplsjM, 



2 M^ 



(6) 



5lf ^ {v,/V2)Tu,J^, 



2 /l/f2 



(7) 



Here a, /3 = 1, 2, 3 (a 7^ (3) denote the quark fiavours. The QFV parameters in the up-type 
squark sector relevant for this study are ^^3^, ^23^-^, 82^^ = (6^^^)* and 823^ = {S^^^)* 
which are the cl — ti, cr — tR, cr — II and cl — tR mixing parameters, respectively. For 
the down- type squark sector we define the QFV parameters as follows (a 7^ /3): 



XdRR 



^ M^^/y^ 



2 /V/f2 



(8) 



XdRL 



{v,/V2)Tr>,jjMl^^Mlpp 



(9) 



The QFV parameters in the down- type squark sector relevant for our study are 5^q^, 5o?^^. 



23 ' "23. _ 

((^3^^)* and ^2^^ = (5^2^)* which are the sl — hi, sr — Sr, sr — hi and sl — iiR 
mixing parameters, respectively. 



5T 



In our analysis we neglect mixing between the first two squark generations due to the 
severe experimental constraints from K meson physics. We also neglect mixing between 
the first and third squark generations focusing on the effects of mixing between the second 
and third generations. We assume all the QFV parameters to be real. These parameters 
are also subject to the experimental constraints given in Table 1. 

Furthermore we impose the vacuum stability conditions for the trilinear coupling ma- 
trices [28] 

\Tuaa\' < 3r^, (Mj„„ + M2„^ + m2), (10) 

\TDaa\' < ^yL{Ml^. + Ml^^ + mi), (11) 

\Tuo.p? < Y^^{M^^^ + M^pp + ml) , (12) 

iTDa^f < Y^^iM^Q^^ + Ml^^ + ml) , (13) 

where a, /3 = 1, 2, 3, a j^ f3; 7 = Max(a, (3) and ml = (m^± + m| sin^ 9w) sin^ /3 — |m|, 
711,2 = ("^H± +"^| sin^ 9w) cos^ (3 — \Tn?z. The Yukawa couplings of the up-type and down- 
type quarks are Yua = V2mu^/v2 = -^ ^^X/j ("^ = u,c,t) and Yoa = V2md^/vi = 
75 m~~cos5 ('^cr = d, s, b), with rriu^ and rrid^ being the running quark masses at the weak 
scale and g being the SU(2) gauge coupling. All soft-SUSY-breaking parameters are also 
assumed to be given at the weak scale. As SM parameters we take mw = 80.4 GeV, 
mz = 91.2 GeV and the on-shell top-quark mass rrit = 173.3 GeV [29]. We have found 
that our results shown in the following are fairly insensitive to the precise value of rrit. 
In addition to the constraints in Table 1 we impose the following limits: 

1. The limit on {rrij^o, tan/3) from the negative search for neutral MSSM Higgs bosons 
decaying into a tau pair (i.e. A^ /H'^ /hP — )■ t'^t") at LHC [30] with rrtjifi being the 
mass of the CP-odd neutral Higgs boson A^ . 

2. The experimental limit on SUSY contributions on the electroweak p parameter [31]: 
Ap (SUSY) < 0.0012. 

3. The LEP limits on the SUSY particle masses [32]: m~± > 103 GeV, m^o > 50 GeV, 
where 'm-± and m.^o are the masses of the lighter chargino and the lightest neutralino, 
respectively. 

We take into account all the constraints listed in this section in all plots presented 
in this article. The constraints on the QFV parameters from B{b — )■ 57), AMb^ and the 
vacuum stability conditions are especially important for this study. For the computation 
of the observables (i.e. physical masses, decay branching ratios, AMb^ and Ap(SUSY)) 
we use the public code SPheno v3.0 [33, 34] . 

3 QFV three-body decays of gluino 

If all squarks are heavier than the gluino and squark generation mixing occurs only be- 
tween the second and third generation, one has the following QFV three-particle decays 





Figure 1: Feynman diagrams for g ^ ctx 



i ■ 



Table 2: Weak scale parameters aX Q = 1 TeV for our prototype QFV scenario, except 
for m^o which is the pole mass (i.e. physical mass ) of A^ . All of Tuaa and Tdqq, are 0. 



Ml M2 M3 fi tan/3 m^o 


139 GeV 264 GeV 800 GeV 1000 GeV 10 800 GeV 












a = /3= 1 


a = /5 = 2 


a = p = 3 






Kc, 


(3150)^ GeV^ 


(3100)^ GeV^ 


(3050)^ GeV^ 






MLp 


(3000)^ GeV^ 


(2200)^ GeV^ 


(2150)^ GeV^ 






ML0 


(3000)^ GeV^ 


(2990)^ GeV^ 


(2980)^ GeV^ 





of gluino into quarks and neutralinos Xi*, i = 1,2, 3, 4, 
9^ shyil, sbx^^. 



(14) 
(15) 



Table 3: Physical masses of the particles in the scenario of Table 2. rriHO is the mass of 
the heavier CP-even neutral Higgs boson H^. 



m^o 


^x°. 


^n 


rris-o 

Xi 




m-+ 
X2 


139 GeV 


281.3 GeV 


1017.9 GeV 


1021.7 GeV 


281.5 GeV 


1022.7 GeV 






rrig 


m^o 


nij^o 


m^o 


mH+ 



975 GeV 121.1 GeV 800.3 GeV 800 GeV 804 GeV 



We will mainly focus on the decays into x5- The corresponding Feynman diagrams for 
the decay g — )■ cix^ are shown in Fig. 1. We have interference between the t and the u 
channel exchanges, as well as between the different Uj exchange diagrams. In particular, 
there can be a strong destructive interference between the ui and Uk contributions, if they 
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Figure 2: (a) Up-type squark and gluino masses and (b) flavour decomposition of Ui (i.e. 
|i?"5p = \cji component!^ (full red line) and l-RigP = \tR component p (dashed blue line)) 
as functions of the QFV parameter 62^^, with the other QFV parameters being zero, for 
the scenario of Table 2. 



are mainly c and t mixtures and their masses are similar. For instance, if ui ~ cos Ota + 
sin 9cji and -Ufc ~ — sin 9iji + cos Ocr then the ui exchange contribution is ~ (T!.^T\ 

whereas the Uk exchange contribution is ~ , T! '^™-, ■ These two contributions almost 
cancel with each other for m,-,, ^ m^ . The suppression of this cancellation requires a large 



"«; 



mass-splitting between the two squarks which can be induced by a large Cf{ — t^ mixing 
term M^^^ (or ^23*^'^) even in case the cr mass parameter M^^^ is similar to the tn mass 
parameter M^^^ . Moreover, in this case one has a very strong cr — Ir mixing. Therefore 
one can expect sizable QFV decay branching ratios for a large mass-splitting and hence for 
large values of ^23^^. For the decays into down-type quarks one has analogous Feynman 
diagrams with the replacements c ^ s, t -^ b, Uj ^ dj. 

There are also gluino three-body decays into charginos, such as ^ -> cb~^, sixf, etc.. 
We will, however, not discuss them explicitly here, although they are included in our 
branching ratio calculations. 

We calculate the three-particle decay branching ratios of the gluino according to the 
diagrams in Fig. 1 and their charge conjugated ones, including the QFV couplings given 
in [6]. As basic SUSY parameters at the weak scale we take Mi, M2, M3, /i, tan/3, m^o, Mq^c 
M^^l^, M^^l^,Tuai3 and Td„/3, which we assume to be real. Here Mi^2,3 are the U(l), 
SU(2) and SU(3) gaugino mass parameters, respectively, and m^o is the pole mass (i.e. 
physical mass) of the Higgs boson A^. We study in detail the QFV scenario based on 
the parameters of Table 2, given at the scale Q = 1 TeV, according to the SPA con- 
vention [35] except for m^o being the pole mass of A^. The scenario of Table 2 satis- 
fies all the constraints listed in Section 2; e.g. for the low energy observables we obtain 




Figure 3: Contours of the QFV decay branching ratio B{g — )■ ctxi) i'^ ^^e S23 - (^23^^ plane 
with the other QFV parameters being zero, for the scenario of Table 2 (solid lines). Also 
shown are the contours of AMb^ = 21.1 ps~^ and 10^xB(6 — )■ 57) = 4.23 (dashed lines). 
The region between the two dashed lines is allowed by all the constraints mentioned in 
Section 2, including those from AMb^ and B(6 — )■ S7). 



AMb, 

S fl~^fl~ 



19.01 ps-S B(6 -> 57) 



5.02 X 10- 



R 



SUSY 

B^TU 



= 3.46 X 10^^ 
0.99, m/jo 



B(6 -^ sl+r) = 1.59 X 
121.1 GeV, Ap (SUSY) = 



10-^ B(6 -^ 
5.70 X 10-5. 



We add to the parameters of Table 2 the QFV parameters 5^3^, 62^^, 51^^^ , '^23'^ ^^ ^^^ ^^ 
52^^i ^23^1 ^tz^ (given also at Q = 1 TeV), and vary them in a range allowed by the con- 
straints listed in Section 2. The physical masses for the case with all the QFV parameters 
being zero are shown in Table 3. They are calculated from the basic MSSM parameters 
at the one-loop level, taking into account the complete flavour structure [33]. We have 
found that these masses are fairly insensitive to the QFV parameters. 

Note that in our case QFV left-right mixing effects, i.e. those due to 82^^, ^23"^) 
^23^^ 5 ^23^^ cannot be significant. We show this for the left-right mixing parameter (Jj-f^. 
Due to the vacuum stability condition (12) we have 17(7321^ < ^(73(^033 + ^^2 
M^33 + Mfj22 ^ 0(10 TeV2), because Yus ^ 1 and mf^ M^^^ ' '>^2 ' 

|T[/32| 



+ ml) 



'C/22- 



fore. 



I KuRL I 
r23 I 



V2 



^ V^^Q33A^&22 



< 



V2_ 
V2 



%33^'tj22 



0.1. Analogously, the parameters 



'^23'^' ^23^^' ^2^^ ^^^ ^ls° constrained to be very small due to the vacuum stability condi- 
tions. Therefore, the most relevant QFV parameters in our study are dJ^i"-, ^2^^, ^23^- 

Fig. 2a shows the physical masses of the up-type squarks Ui, ...,uq as functions of the 
QFV parameter 62^^^, with all the other QFV parameters being zero, for the scenario 
of Table 2. All the constraints mentioned in Section 2 are fulfilled in the shown range. 
Masses of all the down-type squarks di are about 3 TeV in this range. For jf^gl^^l < 0.8 
all squarks are heavier than the gluino. In Fig. 2b we show the flavour decomposition of 
Ml. For |523^^l ^ 0-27 ui is practically a full mixture of cr and t^. For |523^^l ^ 0-88 the 
flavour decomposition of U2 is similar to that of ui with cr and tn interchanged. 



In Fig. 3 we show contours of the branching ratio B(^ -^ ctXi) = B(^ — ?> ctXi) + B(^ — >■ 
ctXi) ill the 62-^ — ^2^^ plane together with contours of 10^xB(6 — )■ S7) = 4.23 and 
AMb^ = 21.1 ps~^, where the other parameters are fixed as in Table 2. The branching 
ratio B(^ — !■ ctx?) can reach values larger than 35%. As can be seen, in the range 
-0.65 < 5^3^ < 0.35, the B(6 -^ 57) and the AMb, constraints are fulfilled. In Ref. [36] it 
is shown that also the p-parameter data may constrain the flavour off-diagonal elements 
of the squark mass matrices, in particular the 5^^ entries. However, in our case the p- 
parameter practically does not give constraints for two reasons. First, we have negligible 
left-right squark mixing implying that the left-squark sector is almost decoupled from the 
right-squark sector. Second, the mass matrices of the left up-type squarks and the left 
down-type squarks are related by the SU{2)l symmetry leading to approximately the 
same masses and mixing matrices for the up-type and down-type squarks. 

In Fig. 4 the branching ratios B(^ — )■ ctxi), B{g — )■ ccxi) ^^^ ^id ~^ t^Xi) ^^^ shown 
as functions of ^23*^'^) "with the other QFV parameters being zero and the other parameters 
fixed as in Table 2. All the constraints mention in Section 2 are fulfilled in the shown 
range. One can see that the QFV decay branching ratio B(^ — >■ ctxi) can reach 40% 
and that in the range 0.6 < |52|^^| < 0.8 the QFV decay branching ratio B{g — > ctx?) is 
even larger than the quark- flavour conserving (QFC) branching ratio B(^ — ?> tixi)- For 
l'^23*^'^l ^ 0-8) the two-body decays into Ui dominate because Ui becomes lighter than the 
gluino (see Fig. 2a). The reason for this large QFV decay branching ratio is as follows: 
For 0.6 < l^g/^^l < 0.8, all squarks other than ui (including down-type squarks) are very 
heavy (see Fig. 2(a)), which leads to the dominance of the Ui exchange contribution in 
the gluino decays. In this 62^^ range the ui, U2 are strong mixtures of c^ and t/j and the 
mass-splitting between ui and U2 is very large, preventing a strong destructive interference 
between the ui and U2 exchange contributions in this ^23*^'^ range (see Fig. 2). This gives 
the large QFV decay branching ratio B{g -^ ctx?)- Note that 'Ui(~ cr + t/j) couples to 
x5(~ -B°) and practically does not couple to X2(~ ^°) ^^^ xf (~ W^) (see Table 2). 
Moreover, X34 and xf are very heavy in the QFV scenario considered here (see Table 3.). 
5° and W^'^ are the U(l) and SU(2) gauginos (the bino and winos), respectively. 

On the other hand, as can be seen in Fig. 3, the dependence of B{g — )■ ctx?) on the 
cl - ti mixing parameter 82^' is much weaker than that on 52^^- This is mainly due to 
the fact that in our scenario the left-squarks {cl^Il) are significantly heavier than the 
right-squarks {cr, Ir) and that the left-squark coupling to Xi (~ bino) is small. Hence, 
the contributions of the left-squark exchanges to B{g ^ ctx?) are suppressed, leading to 
the very small effect of the cl — t^ mixing parameter 62^ on the QFV decay branching 
ratio. As a consequence, for |(523^^| ^ 0.4 the QFV decay branching ratio B(^ —)■ ctxi) is 
smaller than 5% even for larger allowed values of |(5^3^|. 

The gluino can also have QFV decays into down-type quarks with sizeable branching 
ratios if the 62^^ or 5^3^ are unequal to zero. As an example, in Fig. 5 we show a contour 
plot of the QFV decay branching ratio B(^ — )■ sbx^) = B(^ — )■ sbxl) +B{g ^ sbxi) in 
the ^23*^'^ " ^23^ plane for the scenario of Table 2. The dashed contour lines for AMr^ = 
14.5 ps~^ and 21.1 ps~^ show the region —0.38 < 62^^ < 0.12 allowed by all the constraints 
mentioned in Section 2, including the B(6 — t- 57) constraint. Note that in this case the 
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Figure 4: The branching ratios of the decays g -^ cixi + ctxi, g — >■ ccxi and g — )■ ttxi ^^ 
functions of ^23^^ ^^ ^^^ other QFV parameters being zero and the other parameters are 
fixed as in Table 2. o-s 




"23 



cLL 
■'23 



xdRR 
"^23 



Figure 5: Contours of the QFV decay branching ratio B{g -^ sbxi) in the 6. 
plane, with the other QFV parameters being zero for the scenario of Table 2, but with 
the values of M^^^ and M'jj^^ interchanged (solid lines) . Also shown are the contour lines 
for AM^^ = 14.5 ps~^ and AM^^ = 21.1 ps~^ (dashed lines). The region between the 
two dashed lines is allowed by all the constraints mentioned in Section 2, including the 
AMb^ constraint. 



AM^^ constraint is significantly stronger than the B(6 — )■ 57) constraint in the whole 
'^23" ~ '^23^'^ plane. B(^ — )■ shx\) can reach values larger than 30%. The reason for this 
sizable QFV decay branching ratio is similar to that for the large QFV decay branching 
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Figure 6: Contour plot for B{g — )■ ctxi) (solid lines) in the /i — M2 plane for 62^^ = 0.8, 
the other QFV parameters being zero, and the other parameters specified as in Table 2 
with Ml = 139 GeV. Region (A): bino-like LSP region; region (B): wino-like LSP region; 
region (C): higgsino-like LSP region. The point "X" corresponds to our reference scenario 
given in Table 2: M2 = 264 GeV, /i = 1000 GeV. 

ratio B(^ — > ctxi)- The dependence of this QFV decay branching ratio on 5^3^ is again 
much weaker than that on 62^^- 

We also would like to note that we have found a scenario giving QFV three-body decay 
branching ratios B{g — )> ctXi) (or B{g — t- sbxi)) of about 50% for a gluino mass of ~ 1 
TeV, still satisfying all the relevant constraints. In such a scenario, however, the heavier 
squarks have masses of about 6 TeV, while the lightest squark heavier than the gluino 
has a mass of about 1 TeV. 

4 Influence of the neutralino/chargino parameters on 
the QFV three-body gluino decays 

As the squark generation mixing enters into the squark-quark-neutralino/chargino cou- 
plings, here we study how the pattern of the QFV gluino decays depends on the parameters 
of the neutralino-chargino sector. First we show in Fig. 6 a contour plot of the branching 
ratio B{g — )■ ctx?) in the /i — M2 plane for 62^^ = 0.8, the other QFV parameters being 
zero, and the other parameters fixed as in Table 2. In the whole plane rrig ^ 972 GeV. 
As one can see, this branching ratio is larger than 10% for /i > 350 GeV. We indicate the 
regions where the LSP is bino-, wino-, or higgsino-like. The largest QFV decay branching 
ratio B{g — )■ ctx?) is in the bino-like LSP region reaching up to ~ 40%. 

In the following we discuss the dependence of the QFV decay branching ratios of the 
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Figure 7: The /i dependence of the QFV and QFC gluino decay branching ratios for 
(5"^-^ = 0.8, the other QFV parameters being zero, for the scenario given in Table 2. (a) 
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Branching ratios of the decays g — )■ cix^+ctx^, i = 1, ..., 4 as a function of /i. (b) Branching 
ratios of the decays g — ;■ ccxi, g — ?• tixi, g -^ ctx? ' ""^""^ 
function of fj,. 



ctxl and g -^ ctEJ^'' + ctE^'' as a 



gluino on the higgsino mass parameter /i in more detail. In Fig. 7 we show this dependence 
for 62^^ = 0.8, the other QFV parameters being zero, and the other parameters fixed as 
in Table 2. In the whole range nig ^ 972 GeV and the constraints mentioned in Section 2 



cO 



1, 



are satisfied. The branching ratios of the QFV gluino decays g — )■ ctXi + ctx 
are shown in Fig. 7a. For fi > Mi the bino component of x? is increasing and hence 
the branching ratio of ^ — )■ ci{ct)xi increases. For /i > 700 GeV this branching ratio is 
about 36%, being roughly a factor of 2 larger than that for g — > ttXi- For l/^l ^ ^i; 
Xi and X2 ^^6 essentially higgsinos, and the branching ratio of ^ — ;■ ci{ct)xi is less than 
10%. The reason is that the QFV decays into the higher neutralinos and the charginos 
become more important. Note that |/i| < 100 GeV is excluded by the LEP limit on 

can contain contributions from 



the m-± 



Furthermore, the final state ctE'™** 



+ ctEJ^'' 



the higher neutralino modes ct(ct)Xi , "^ > 2, with the invisible decays of Xi ~^ Xi'^^^ 
where i?™** is missing transverse energy. Therefore, we show in Fig. 7b a plot where these 
contributions are included (dashed red line). One can clearly see that for /j, < 320 GeV the 
contributions of the invisible decays of the higher neutralinos (see Fig. 7a) are important. 
For comparison in Fig. 7b the QFC branching ratios B{g — )■ ccXi) and B{g — )■ tixi) are 
also shown. 

The situation is quite different if Mi > M2. In this case, for /i > M2 the Xi is 
essentially a wino which does not couple to cr and Ir leading to a small branching ratio 
of ^ — ;■ ci{ct)x\. On the other hand, a large branching ratio of ^ — )■ ci{ct)x^ can be 
expected, because for |/i| > Mi, X2 becomes bino-like. 

Next we discuss the dependence of the QFV gluino decay branching ratios on the 
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SU(2) gaugino mass parameter M2. In Fig. 8 we show this dependence fixing Mi = 
264 GeV, fi = 600 GeV, 62-^^ = 0.8, the other QFV parameters being zero, and the 
other parameters fixed as in Table 2. In the shown range rrig ^ 972 GeV and all the 
constraints mentioned in Section 2 are satisfied. As just explained above, in Fig. 8a, for 
M2 < Ml the branching ratio of the decay g — )■ ci{ct)xi is almost zero, while that for 
g — !■ ci{ct)x2 is large (~ 24%). For M2 > Mi the roles of x? and X2 ^^^ interchanged 
and therefore the decay g — ?■ ci{ct)xi becomes dominant with a branching ratio of about 
25%. Note that the range M2 < 100 GeV is excluded by the LEP chargino mass limit 
mentioned in Section 2: m-± > 103 GeV. At M2 ~ u there is again a level crossing. The 
X4 is higgsino-like for M2 < fi and becomes wino-like for M2 > fi. 

In Fig. 8b the branching ratios for the decays g — ?■ cb{cb)x^ as a function of M2 are 
shown. The level crossing of Xi and X2 at M2 ~ /i is clearly seen. 

In the following we discuss in more detail a typical scenario with a higgsino-like LSP 
(/i < Mi,M2) and one with a wino-like LSP (M2 < Mi,/i). As an example for the 
higgsino-like LSP scenario, we choose the parameters as given in Table 4(a), with the 
squark mass parameters as in Table 2. We fix the QFV parameter 62^^ = 0.8 and 
the other QFV parameters equal zero. In this scenario all experimental and theoretical 
constraints mentioned in Section 2 are fulfilled. The relevant masses for the neutralinos 
and the charginos are given in Table 4(b). We show the most important QFV decay 
branching ratios in Table 4(c). In this scenario X12 and Xi are almost higgsinos and 
hence their couplings to Mi(~ cr + tR) are significantly enhanced by the large top-quark 
Yukawa coupling, which results in the sizable branching ratios of the QFV decays into 
Xi 2 and Xi . Since B(x2 -^ Xi^^) is relatively large (= 18.4%), the sum of the branching 
ratios for the g decays into the final states (ct + i?™**) and (ct + i^^™**) is sizable (= 9.0%). 
The leptonic xf decays xf -^ Z^^^'/iXi and xf -^ e^VeXi have a branching ratio of 
approximately 11.2% each, because W exchange dominates. Therefore, from the gluino 
decays into Xi 0116 gets final states hc{hc) -|-/i^(e^) + E1^^^ with a branching ratio of 5.0%. 
This has to be compared with the expectation from MFV which is of the order of 10~^, 
as it is proportional to \Vch\^. 

Next we discuss a scenario where the LSP is wino-like, with the parameters as given 
in Table 5(a), the squark mass parameters as in Table 2, 82^^ = 0.8 and the other 
QFV parameters being zero. Again, all the constraints are satisfied. The masses of 
the neutralinos and charginos are given in Table 5(b). The relevant QFV gluino decay 
branching ratios are shown in Table 5(c). As in this case B(x2 -^ x\^^) = 4.1%, the sum 
of the branchung ratios for the final states {ci+ E^'^^) and (ct -|- i^^™"*) is 6%. As B(x^ — >■ 
fJ'^^f.Xi) = Hxt -^ e+PeXi) = 13.2%, one has B(^ -^ bcfi+u^Xi) = B(^ -^ k/i-i/^x?) 
= B{g — !■ bce'^h'eXi) = B{g — )■ bce'UeXi) = 0.8%. Hence, the signature be (or be) plus a 
lepton (/i^*" or e^) plus E!^^^ has a probability of about 3%. 

Summarizing the discussion of this section we can say that the branching ratios of the 
QFV three-particle gluino decays depend not only on the generation mixing in the squark 
sector, but also quite strongly on the parameters of the neutralino/chargino sector. 
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Figure 8: The M2 dependence of the QFV gluino decay branching ratios for Mi = 
264 GeV, /i = 600 GeV, 62^^ = 0.8, the other QFV parameters being zero, and the 
other parameters fixed as in Table 2. (a) Branching ratios of the decays g — )■ cix^ + ctx^, 
i = 1, ...,4 as a function of M2. (b) Branching ratios of the decays g — )■ d>x^ + cbxt, 
k = 1,2 as a function of M2. 



Table 4: Weak scale parameters at Q = 1 TeV (except for ?n^o being the pole mass), the 
corresponding neutralino and chargino masses and some important branching ratios for a 
scenario with a higgsino-like LSP, where rrig = 972 GeV. 

(a) Weak scale parameters at Q = 1 TeV. 



Ml M2 M3 fi tan/3 itlj^q 


139 GeV 264 GeV 800 GeV 120 GeV 10 800 GeV 





(b) 


Neutralino and chargino masses. 




Xi 


"^x§ 


X3 


X4 


"^xt 


m-+ 
X2 


87.0 GeV 


133.1 GeV 


158.3 GeV 


310.1 GeV 


109.3 GeV 


310.1 GeV 





(c) Important branching ratios. 




B(5 ^ cixl) 


B{~g^cix'2) B{~g^bcxt) B(x^ ^ xM 


B{xt ^ f^+i^f.xl) 


3.4% 


6.1 % 11.2 % 18.4 % 


11.2 % 
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Table 5: Weak scale parameters a.t Q = 1 TeV (except for tti^o being the pole mass), the 
corresponding neutralino and chargino masses and some important branching ratios for a 
scenario with a wino-like LSP, where tti^ = 972 GeV. 





(a) Weak scale parameters at Q = 1 TeV. 




Ml 


M2 M3 /i tan/3 


m^o 


400 GeV 


300 GeV 800 GeV 350 GeV 10 


800 GeV 



(b) Neutralino and chargino masses. 



Xi 


771^0 
A2 


X3 




m-+ 

Xi 


X9, 


275.5 GeV 


362.6 GeV 


376.4 GeV 


433.1 GeV 


280.0 GeV 


407.1 GeV 





(c) Important branching ratios. 




B(5 ^ ctx?) 


BQ^cix^) B(5^6cxr) B(x^^xM 


B(xr ^ Z^+'^mX',') 


2.5% 


6.3 % 5.8 % 4.1 % 


13.2 % 



5 Measurability of the QFV gluino three-body de- 
cays 

We calculate the relevant gluino production cross sections at leading order using the 
WHIZARD/O'MEGA packages [37, 38] where we have implemented the model described 
in Section 2 with squark generation mixing in its most general form. We use the CTEQ6L 
global parton density fit [39] for the parton distribution functions and take Q = m^^ +mp2 
for the factorization scale, where pi and p2 are the sparticle pair produced. The QCD 
coupling as{Q) is also evaluated (at the two-loop level) at this scale Q. Due to the heavy 
squarks in our reference scenario of Table 2, the dominant gluino production process at 
LHC is pp — )■ ggX, where X contains beam-jets only. For the scenario of Table 2 the 
corresponding cross section is practically independent of 62^^ and is about 170 fb (3 fb) 
at v^ = 14 TeV (7 TeV). (Note that m~g = 975 GeV and 972 GeV for 6^-^^ = and 
0.8, respectively.) The sum of the cross sections of the other gluino production processes, 
such as pp — >■ gqiX, gq2X, gXiX, gx2^, is two orders of magnitude smaller than that of 
pp -^ ggX. 

In Fig. 9 we show the signal rates due to pp — )■ ggX, with X containing beam-jets only, 
at y/s = 14 TeV, where at least one of the pair-produced gluinos decays as ^ — )■ ci{ct)x1, 
as a function of 62^^ for the scenario of Table 2. All the constraints mentioned in Section 2 
are satisfied in the range |523^^| < 0.85. The rate of the final state ccttE^^'^ X , produced 
in the case when both gluinos decay like g — > cix\-, reaches 7 fb for I'Jgl^'^l = 0.8 (full blue 
line), yielding 700 events for an integrated luminosity of 100 fb~^. The charge conjugated 
final state ccttE^^^X has the same rate. The full red line shows the rate for the QFV 
case where one gluino decays as ^ — )■ cix\ and the other one as g ^ ctxi plus the QFC 
case with one g — )■ ccxi and the other g — >■ tixi- The dashed red line presents the rate 
for the QFC case only. One can see that the QFV signal has a rate of about 14 fb for 
0.8. The green lines show the case of one gluino decaying as ^ — )■ cix^, and the 
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Figure 9: Signal rates for pp — )■ ggX at ^/s = 14 TeV where at least one of the gluinos 
decays as ^ — ?■ ct(ct)xi, as a function of ^23^'^ with the other QFV parameters being zero 
and the other parameters fixed as in Table 2. Shown are the rates for the final states 
with cciiE^'' (full blue hue), cctiE^'' (QFV + QFC) (full red line), cctiE^'' (QFC only) 
(dashed red line), ctttE^"^^ (dashed green line), ccciE^"^^ (full green line). 



other one as ^ — )■ ccxi (full green line) or g ^ tixi (dashed green line). These rates reach 
25 fb and 10 fb, respectively, for |52|^^| = 0.8. The charge conjugated final states have 
the same rates. 

A characteristic feature of a three-particle decay of the gluino is the invariant mass 
distribution of the two produced quarks. We calculate the invariant mass distributions 
dr(^ -^ UjUkXi)/{Ttot{9) dMu^Uk 






system UjUk, 



(Pu, 



where Muu^. is the invariant mass of the two-quark 



+ Pu,. 
:uRR _ 



In Fig. 10 we show these distributions for g — )■ 
ttxi, ccXi, ctxi + ctxi with 62^^ = 0.8 and the other QFV parameters being zero for 
the scenario of Table 2. In contrast to the case where the squarks are lighter than the 
gluino and the gluino decays via real squarks [11], no edge structure appears. However, 
the thresholds and the shapes of the distributions are very different. The endpoint is at 
[nig — m^o). The thresholds are at 2mc,mc + rrit and 2mt, respectively. Measuring these 
distributions could be helpful for separating the QFV decays into cixi + ^tx? from the 
QFC decays. 

A typical gg production event has at least four laige-px jets and large E'™*. An event 
with a QFV gluino decay g — )■ ctxi(ctxi) should contain at least one top (anti-top) quark 
in the final state, which must be identified. This is possible by using the decay t — )■ bW^ 
with the W^ decaying into two jets. For this purpose a special method was proposed 
in [40] . Charm tagging would be extremely helpful for a clear identification of the QFV 
gluino decay g — )■ ctxi{cixi)- If this is not possible one could search for the decay g — >■ 
QtXiiQtXi) y Q ¥" t- Typical signatures of QFV g pair events are: t(or t) +3 jets + -E'T*'* + ^) 
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Figure 10: Invariant mass distributions of two up-type quarks from the decay g — )■ UjUkXi, 
with 62^^ = 0.8, the other QFV parameters being zero, and the other parameters fixed 
as in Table 2. 



t + t{oTt + i)+2 jets + E'^^' + X and t + 1 + 1 (or t + 1 + 1) + 1 jet + E^^' + X , where X 
contains beam-jets only. Note, that the signal events t + t (or i +t) + 2jets + E'™* -|- X 
can practically not be produced in the MSSM (nor in the SM) with QFC 

The rate of a possible SUSY background from pair production of squarks, such as 
pp — )■ qqX with q — )■ ex?, q — )■ ixi is much smaller than that of the signal of gluino pair 
production due to the larger squark masses. As shown in the SUSY searches by ATLAS 
and CMS [12, 13], the SM backgrounds, such as QCD multijets, W^ + jets, Z^ + jets, ti 
and single top production, can be strongly reduced by appropriate selection cuts. 

6 Conclusions 

We have studied QFV decays of gluino within the MSSM in the case that all squarks are 
heavier than the gluino and, hence, the gluino has only three-particle decays. Starting 
from the most general squark mass matrix, we have assumed mixing between the second 
and the third squark generations in the up and down sectors. We have taken into account 
all relevant experimental constraints from SUSY particles and Higgs searches as well as 
from precision data in the B meson sector. Furthermore we have respected the vacuum 
stability conditions for the trilinear coupling matrices. It has turned out that of all 
QFV parameters the parameters 62^^, ^23^^ ^^^ ^23 V^^Y the most important role in our 
study. We have concentrated on the QFV decays g — > ct(ct)xi and g — )■ hs{hs)x\ which 
presumably have the clearest signatures for the presence of QFV in the MSSM. We have 
studied these within a prototype scenario with gluino mass rrig ^ 1 TeV and squarks with 
masses between ~ 1 TeV and ~ 3 TeV, where the lightest squarks are c/j — tfi mixtures 
(or SR — bji mixtures). These QFV decay branching ratios can reach up to 40% (35%). We 
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have paid special attention to the dependence of the QFV decay branching ratios on the 
chargino/neutrahno parameters. In this context we have considered three cases where the 
hghtest neutrahno is bino-, wino- and higgsino-hke, respectively. We have found that the 
QFV decay branching ratios depend strongly also on the chargino/neutralino parameters. 
As in our scenario the squarks are heavier than the gluino, the dominant gluino pro- 
duction process at LHC is gluino pair production, pp -^ ggX. We have calculated the 
rates for the various signatures stemming from QFV gluino decays as well as the invariant 
mass distributions of the two-quark system ci+ ct, cc, ti in the final states. We have found 
that the rates of the resulting QFV signatures, like pp — )■ tcccE'!^'^^ X and pp — > ttccEl^'^^ X , 
are significant at LHC. This could have an important influence on the search for gluinos 
and the determination of the basic MSSM parameters at LHC. 
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